Universal nature of particle displacements close to glass and jamming transitions 
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We examine the structure of the distribution of single particle displacements (van-Hove function) 
in a broad class of materials close to glass and jamming transitions. In a wide time window com- 
prising structural relaxation, van-Hove functions reflect the coexistence of slow and fast particles 
(dynamic heterogeneity). The tails of the distributions exhibit exponential, rather than Gaussian, 
decay. We argue that this behavior is universal in glassy materials and should be considered the 
analog, in space, of the stretched exponential decay of time correlation functions. We introduce 
a dynamical model that describes quantitatively numerical and experimental data in supercooled 
liquids, colloidal hard spheres and granular materials. The tails of the distributions directly explain 
the decoupling between translational diffusion and structural relaxation observed in glassy materials. 

PACS numbers: 02.70.Ns, 64.70.Pf, 05.20.Jj 



The slow dynamies of disordered materials close to 
glass and jamming transitions is characterized by just 
a few universal features [H : dramatic dynamical changes 
upon mild changes of control parameters (temperature, 
density), broad distribution of relaxation times leading to 
stretched exponential decay of time correlation functions, 
and spatially heterogeneous dynamics [3] . Here we argue 
that the detailed structure of the distribution of parti- 
cles displacements (van-Hove function Q ) constitutes an 
additional universal signature of glassy dynamics. We 
show that, for timescales corresponding to structural re- 
laxation, the self part of the van-Hove function has fat 
tails that are well-described by an exponential, rather 
than a Gaussian, decay. We provide a broad range of nu- 
merical and experimental data, physical arguments, and 
a dynamical model to support this claim. 

The non-Fickian character of single particle displace- 
ments in materials with glassy dynamics is well-known 0, 
H, 0, 0, H, S [13, EHj Ell ■ correlation functions decay 
non-exponentially, mean-squared displacements exhibit 
at intermediate timescales a subdiffusive plateau, van- 
Hove distributions are non-Gaussian. This affects trans- 
port properties since translational diffusion is decoupled 
from structural relaxation , leading to an anomalous 
relation between timescales and lengthscales [ij] . Virtu- 



ally all glass theories address the stretched decay of cor- 
relation functions, but comparatively much less attention 
has been paid to the detailed shape of the self-part of the 
van-Hove function [ll|, [H, E [IE [S [Hi , although new 
techniques now directly access this quantity in different 
materials [3, [1, [1, [3- non-Gaussian, "heterogeneous" 



character is often discussed in qualitative terms [16|, ll7| , 
and quantitative measures focus on the distribution kur- 
tosis (non-Gaussian parameter [^) which contains, how- 
ever, very indirect information about its shape. Devia- 
tions from Gaussian behavior are usually ascribed to dy- 
namic heterogeneity [1|, i.e. to the presence of particles 
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FIG. 1: Time evolution of the self part of van-Hove func- 
tions for silicon atoms in silica, Lennard- Jones particles, hard- 
sphere colloids and grains (open circles), fitted with the 
model in Eq. ([2]) (full lines). They exhibit a Gaussian cen- 
tral part and a fat, exponential tail, (a) T = 3000 K and 
t G [27, 1650] ps. (b) T = 0.435 and t € [75 ■ 10^,41 • 10®]. (c) 
= 0.517 and t £ [90, 1008] s. (d) ip = 0.84 and t e [10, 1000] 
cycles, (a) and (b) show the distributions of |r(t) — r(0)|, (c) 
and (d) the distributions of x{t) — x{0). 



that are substantially faster or slower than the average. 
We argue that van-Hove functions contain quantitatively 
relevant information about the relaxation of glassy mate- 
rials, and that its functional form is simple and universal, 
just as the stretched exponential decay of time correla- 
tion functions. Glass theories should therefore treat both 
phenomena on an equal footing. 

We present our central observations in Figs. [H [5] which 
show the self part of the van-Hove functions for a sil- 
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FIG. 2: Temperature evolution of the self part of the van- 
Hove function for t — Tc for oxygen atoms in silica at T = 
3580, 3200, 3000 and 2715 K (from left to right) and Lennard- 
Jones particles at T = 0.5, 0.47, 0.45 and 0.435 (from left 
to right). The exponential tail gets more pronounced at low 
temperatures. This trend is smaller in silica than for Lennard- 
Jones, as is the translational decoupling. 



ica melt [lOj, a binary Lennard- Jones (LJ) mixture [20|, 
a dense suspension of colloidal hard spheres and a 
slowly driven dense granular assembly ^. It reads Q 



Gs{r,t) = {6{r-\r,{t)-rm)), 



(1) 



where ri{t) denotes the position of a particle i (molecule, 
colloid, grain) at time the brackets indicate an en- 
semble average. For technical reasons, experiments 
sometimes record the one-dimensional version of p]). 
Gs{x,t) = {S{x — \xi{t) — a;i(0)|)), where Xi{t) is the 
projection of ri{t) on a given unit vector. For our pur- 
poses, the difference between the two functions is irrele- 
vant. For all systems we find that Gs{r,t) has the same 
structure over a broad time window comprising the struc- 
tural relaxation. Most of its statistical weight is car- 
ried by particles that have barely moved, r < cr, but 
a "fat" tail extends to much larger distances, r > <t, 
where a is the particle diameter. The small r behavior 
is not far from a Gaussian distribution, corresponding to 
quasi-harmonic vibrations, but the large distance decay 
appears linear in Fig. [1] i.e. Gs oc exp(— r/A(t)). For a 
Fickian particle, one expects instead a Gaussian decay, 
Gs oc exp{—r^/{4Dst)), where Dg is the self-diffusion 
constant. Although the existence of a fat tail in GJr, t) 
was recognized before 0, i, S, 0, i, [H, E [H, [13, El , its 
non-trivial functional form and universality went largely 
unnoticed. 

The exponential tail extends to larger distances when t 
increases, but \{t) grows very slowly with t. Eventually, 
at very large times, a crossover to Fickian behavior is ob- 
served, see the latest time in Fig.[l]b. For the LJ system 
the crossover takes place for t « 30tq [21| , where Tq (T) is 
defined from the time decay of the self-intermediate scat- 
tering function, Fs (q, t) , the Fourier transform of Gs {t, i) . 
Our observations correspond to times that are shorter 
than this crossover. In Fig. [21 we present the evolution 
of Gs {t, t) when the glass transition is approached, keep- 
ing t fixed to Tq (T) . Clearly the shape of Gs {r, t) remains 




FIG. 3: Four trajectories of duration Tq for particles con- 
tributing to the center of the distribution ("slow") and to 
the tail ("tail") taken from simulations of the LJ liquid at 
T = 0.435 ("LJ") and for random walkers described by ([2]). 



unchanged, but the tail becomes more pronounced closer 
to the glass transition. 

These observations confirm that van-Hove functions 
can be qualitatively described as the superposition of two 
families of particles: localized particles contributing to 
the Gaussian central part and mobile particles contribut- 
ing to the exponential tail. Evidence has recently been 
given that the distinction between mobile and immobile 
particles cannot be explained on a structural basis [2^ . 
We therefore seek a purely dynamical explanation. In 
Fig. [3] we present representative trajectories of duration 
Ta for mobile and immobile particles in the LJ system. 
Similar pictures have been presented before 0, Q- A 
large fraction of the particles perform localized, vibra- 
tional motion around their initial positions, as in a disor- 
dered solid. These "slow" particles are "caged" by their 
neighbors and contribute to the quasi-Gaussian central 
part of Gs {r, t) . 

More interesting is the behavior of the particles con- 
tributing to the tail. On top of the localized vibra- 
tions, these particles perform a (distributed) number of 
quasi-instantaneous "jumps" . This suggests that par- 
ticles perform a form of continuous time random walk 
(CTRW (H) [nl, [H, [H m. From direct inspection 
of the trajectories we note that the size of the jumps is 
distributed, and represents on average only a small frac- 
tion of the particle size, implying that jumps probably 
result from cooperative events involving a large number 
of particles moving by a small amount [2^. Regarding 
timescales, trajectories also reveal that the time of the 
first jump after t = is distributed. This observation di- 
rectly implies that van-Hove functions can be described, 
for t > 0, as a superposition between particles that have 
jumped, and those that have not 
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[151 . We insist 
that this coexistence is dynamically generated and we 
will avoid the assumption of a material being composed 
of two dynamically distinct families of particles [3, [2^ . 

The final empirical observation from Fig.[3]is that once 
a particle has managed to make a jump, it very likely 
makes one or several additional jumps during the rest 
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of our observation time. We believe that this results 
from dynamic heterogeneity. Spatial clustering of parti- 
cles with correlated dynamics implies indeed that it takes 
a very long time for a particle belonging to a slow region 
to become mobile. But when this happens, the particle 
then likely belongs to a mobile region, which enhances 
considerably its probability to move further. Different 
timescales for initial and subsequent moves is an exact 
result, for the reason mentioned above, in kinetically con- 
11, [23. 



strained models 



It is likely a generic consequence 
of the presence of spatially heterogeneous dynamics. 

These features endow particle diffusion in glassy mate- 
rials with specific properties. We now introduce a model 
which incorporates these empirical observations with as 
few free parameters as possible. We significantly extend 
the work of Ref . [11] , which analyzed self-diffusion in ki- 
netically constrained models, to describe off-lattice re- 
alistic models and experiments. The system is viewed 
as an assembly of dynamically indistinguishable parti- 
cles, compatible with structural homogeneity. We as- 
sume solid behavior at short times. In between jumps, 
particles perform on a microscopic timcscalc a Gaus- 
sian exploration of their environment with the distri- 
bution Uihir) = (27r^2)-3/2exp(-rV2^2). Wc assume 
that particles perform jumps with a size sampled from 
/j„^p(r) = (27rd2)-3/2exp(-r2/2d^). We similarly as- 
sume simple forms for time distributions. The time of 
the first jump is drawn from 4>i{t) = t^^ exp(— i/ri). We 
then assume that subsequent jumps arise with higher fre- 
quency, using the distribution 4'2{t) = expi—t/r-p), 
with T2 < Ti. It is now a simple task [2J| to express 
the van-Hove function, Gs{r,t) = ^'^^Qp{n,t)f{n,r), 
where p{n, t) is the probability to make n jumps in a 
time i, and f(n, r) the probability to move a distance r 
in n jumps. These probabilities involve convolutions and 
are more easily expressed in the Fouricr-Laplacc domain, 
(r, t) (g, s). We obtain 

Gs{q, s) = /vib(g)<i>i(^s) + /(g)/vib(g) , ^'^f'ft'} , : (2) 

1 - 02(s)/(g) 

where (/'i,2(s) = 1 - s$i,2(s) and f{q) = /vib(g)/jump('7)- 
The Montroll- Weiss equation [2J] is recovered when (pi = 
4>2 and vibrations are not considered, fvibil) = 1- The 
result in Eq. ^ is valid for any choice of distributions 
(/vib, /jump, 4'2) ■ Hcrc, WC restrict to simple choices 
(Gaussian and exponential distributions) to emphasize 
the universality and physical origin of our results and to 
introduce as few free fitting parameters as possible: (£, 
d, Ti, T2). Equation ([2]) makes very transparent the fact 
that Gs is the superposition of localized particles, and 
mobile particles. We show below that the second term 
produces a tail that is close to exponential and arises from 
particles which have performed one or several jumps. 

For the four systems considered in Fig. [T] we have used 
Eq. ^ to fit the self part of the van-Hove functions, as 
shown with full lines. The fits evidently match the data 
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FIG. 4: Comparison of the amount of translational decou- 
pling in silica and LJ systems measured in MD simulations 
through Rdcc (open symbols) , and of the ratio a — ti/t2 (filled 
symbols) employed to fit the data with the model The 
agreement is excellent. 



very well. In practice, we sought the set of parameters 
that allows for data fitting on the largest time window 
comprising structural relaxation. We find that fitting 
several times fixes the set of parameters with little ambi- 
guity, while multiple choices remain possible when fitting 
data for a single t. In Fig. [3] we present numerically gen- 
erated trajectories of the generalized CTRW model of 
using the parameters used to fit the LJ data, leaving 
us with no doubt that such a model captures the main 
qualitative aspects of the real trajectories. As expected, 
we find that both cage and jump sizes represent only a 
fraction of the particle size, and are very weakly depen- 
dent on the control parameters. As a rule of thumb we 
find d « 2£. For instance {d/a, i/a) is (0.1, 0.051) in 
coUoids, (0.15, 0.06) in grains, and (0.35, 0.15) in the 
Lennard- Jones. Moreover, our results for d^ agree well 
with plateau values directly measured in mean-squared 
displacements for all systems. 

These results imply that the changes observed in Fig. [2] 
are mostly due to a change in the balance between ti and 
T2. We naturally find that, in a first approximation, both 
times track the structural relaxation of the system. More 
interesting is the evolution of their ratio a ~ ti/t2, re- 
ported in Fig. [31 Unfortunately, we had not enough data 
for grains and colloids to report accurate estimates of a in 
these systems and thus we concentrate on the two numer- 
ical models. In order to account for the increasingly fat 
tails reported in Fig. [21 a has to grow significantly when 
temperature decreases, as expected from the above dis- 
cussion. The grow th of a directly impacts on transport 
properties [l5, 27|. Within the model ^ the time decay 
of Fs{q, t) at large q (small distance) is governed by 0i(t), 
so that Tq, ~ Ti . Fickian diffusion is recovered when the 
average number of jumps becomes large. It is easy to 
show from the g, s ^ hmit of © that Dg ^ {f'+d^)/T2, 
so that the product between self-diffusion constant and 
structural relaxation time scales as DsTn ~ a. Our model 
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therefore makes a direct prediction about translational 
decoupling. We have measured the normaUzed prod- 
net, i?dcc(r) = D,{T)T^iT)/D,iTo)T^{To) (where Tq is 
a high temperature) for Lennard-Jones, siUcon and oxy- 
gen atoms directly in numerical simulations, see Fig. |31 
This ratio is i?doc = 1 at high temperatures, and becomes 
i?dec > 1 whenever translational decoupling occurs 14 1. 
For three types of particles with different degrees of de- 
coupling, we find quantitative agreement between a ob- 
tained from fits of the self part of the van-Hove function 
and the decoupling i?dec directly measured in the simu- 
lations. Thus, Fig.|4]gives strong support to our physical 
description and empirical modeling of self-diffusion close 
to a glass transition, and provides a quantitative link be- 
tween dynamic heterogeneity and decoupling. 

Why are the tails of the distributions described by an 
exponential decay? Non-Gaussian decay is in fact present 
in the original CTRW model when distances outside the 
realm of central limit theorem are considered. These tails 
are enhanced, and hence more easily detectable, when 
a > 1 and decoupling occurs. Consider the case a = 1, 
£ = in Eq. Inverting the Laplace transform yields 
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dg[e*-^(«) - l]gsin(gr), (3) 



where Go{r, t) = (5(r)$i(t) and i = t/ri. We then expand 
the exponential in ([3]), integrate each term and convert 
the sum into an integral to get 



Gsir,t) = Goir,t) + 



4d3 



dn- 



=-/(») 



(4) 



with /(n) = n\nn — n\nt — n + r^/(8o?^n). The large 
distance limit of (|4]) is evaluated using a saddle point 
approximation, 



Gsir,t) 



3/2g-t 



,-r[y-l/y']/2d 



(rd)3/2yrTY2 



(5) 



where Y satisfies Y^ expY^ = r^/{2dt)^; Y^ 2\og{^^) 
for large r. Thus, we find that Gs{r,t) decays exponen- 
tially (with logarithmic corrections) at large r. Inter- 
estingly this expansion can be obtained independently 
of the actual shape of the distributions, establishing its 
universality. Considering that the tail of Gs stems from 
particles that have performed a number of jumps larger 
than average, one finds p{n,t) ~ it/nY and f{n,r) ^ 
g-r /(8d n) ^ yielding an expression similar to 

We have reported the existence of a new universal fea- 
ture characterizing the dynamics of materials close to 
glass and jamming transitions, seen in the structure of 
the distribution of single particle displacements which ex- 
hibits exponential decay at large distances. We argued it 
is a generic consequence of the existence of spatially het- 
erogeneous dynamics, which has profound consequences 
on transport properties. Our results apply to a wide va- 
riety of materials from atomic liquids [l| to biophysical 



materials @ and grains . We strongly encourage more 
systematic experimental exploration of particle displace- 
ments in amorphous materials with slow dynamics. 
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